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ABSTRACT 

We determine the radial abundance gradient of Cl in the Milky Way from Hii 
regions spectra. For the first time, the Cl/H ratios are computed by simply adding ionic 
abundances and not using an ionization correction factor (ICF). We use a collection of 
published very deep spectra of Galactic H ii regions. We have re-calculated the physical 
conditions, ionic and total abundances of Cl and O using the same methodology and 
updated atomic data for all the objects. We find that the slopes of the radial gradients 
of Cl and O are identical within the uncertainties: —0.043 dex kpc“^. This is consistent 
with a lockstep evolution of both elements. We obtain that the mean value of the 
Cl/0 ratio across the Galactic disc is log(Gl/0) = —3.42 ± 0.06. We compare our 
Cl/H ratios with those determined from CP+ abundances and using some available 
ICF schemes of the literature. We find that our total Cl abundances are always lower 
than the values determined using ICFs, indicating that those correction schemes 
systematically overestimate the contribution of C1+ and Cl^’*' species to the total Cl 
abundance. Finally, we propose an empirical ICF{Cl'^~^) to estimate the Cl/H ratio 
in Hii regions. 

Key words: ISM: abundances - H ii regions - Galaxy: abundances - Galaxy: evolu¬ 
tion 


1 INTRODUCTION 

Chlorine (Cl) is located in Group VIIA of the periodic ta¬ 
ble. The elements of this group are known as halogens, 
which also includes fluorine, bromine and iodine. Cl is be¬ 
tween the 17t h and 20th most ab undant element in the 
Solar System llAsplund et al.l [20091) and its most common 
isotopes are ®®C1 and ®^C1. These two isotopes are pro¬ 
duced durin g both hydrostatic and explosive oxy gen burn¬ 
ing in stars llWooslev fc Weav^l 19951 : iPagelll 19971 ). although 
some contribution may be produced by Type la supernovae 
llNomoto et al.l 19971 ). IClavtonI (l2003l ) indicates that once ®‘*S 
is formed during oxygen burning, a single additional proton 
capture yields ^®C1. In the case of ®^C1, most of it is pro¬ 
duced by radioactive decay of ®^Ar, which is formed by a 
single neutron capture by ®®Ar. These proce sses seem to ac - 
count for the solar ®®C1 and ®^C1 abundances llClavtonl2003l) . 
Since lines of Cl are very rare in stellar spectra, t here are 
very few determinations of its abundance in stars llClavtorJ 
I2OO3I ). In addition, the Cl abundance is difficult to calcu- 
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late in diffuse and dense interstellar clouds lljenkinsll2009l ). 
Therefore, most of the information we have about Cl abun¬ 
dances comes from the analysis of emission lines in ionized 
nebulae. 

The determination of Cl abundances in H ll regions and 
planetary nebulae (PNe) is normally based on the measure¬ 
ment of the intensity of [Cl ill] doublet at 5517 and 5537 A, 
which are by far the brighest lines of Cl in optical spectra. 
For the common ionization conditions of Hll regions, we ex¬ 
pect to have only significant amounts of CH, CF^ and CF"*". 
Given the similarity between the ionization potential (IP) 
of He+ and CF+, we do not expect a significant amount of 
CF'*' in H II regions with undetected He ll emission lines. The 
brightest lines of CH and CF^ are [Clll] 9123 A and [Cliv] 
8046 A, respectively and, for example, in the Orion Neb¬ 
ula are only between 4% and 8% the intensity of the [Cl ill] 
lines. Therefore, the detection of [Cl ll] and [Cl iv] lines may 
only be possible in deep spectra of bright objects that also 
include the far red (until about 9100 A). Moreover, the red¬ 
dest part of the optical spectrum is heavily contaminated 
by intense sky lines so the use of high spectral resolution 
or an exquisite sky emission removal is mandatory in order 
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to separate these faint nebular lines. In the usual situation 
of only having [Cl ill] lines in the spectra we must adopt an 
ionization correction factor (ICF) to estimate the total Cl 
abundance. In this case, the ICF{Cl^'^) is a multiplicative 
factor to transform Cl^’^/H'^ ratios into Cl/H ones. It de¬ 
pends on the ionization degree of the nebula, generally on its 
/O^'^ ratio or a combination of the 0^/0^'^ and 
ones. For H ll region studies, there are different expressions 
for the JCF (Cl^~*~) available in the li terature. For exam¬ 
ple, [PeimberF&^^^rres^Peimbe^ lll977ll proposed a relation 
based on the similarity of the ioniz ation potential of ionic 
sp ecies of Cl, O and S but others, as iMathis fc Ro^ (1 199 il l 
or lizotov et al.l ll2006ll . obtained fitting expressions based on 
photoionization models. A very co mplete ICF scheme for 
PNe h as been recently presented bv iDelgado-Inglada et aP 
ll2014h . That work includes an JCT( Cl^'*~/0^) that can be 
appli cable for H ll region conditions (IPelgado-Inglada et al.l 

[2^. 


The radial abundance gradient of Cl in the Milky 
Way has only been studied in papers devoted to PNe 
but never in works based on H ll region observa¬ 
tions. I^undez-Abaiisfc Mn^ll lll98lt) . lMaciel fc Chiappinil 
lll994h and Henry et al. 1 2004l have determined the radial 
gradient of Cl in the Milky Way from abundance data for 
large samples of PNe, finding slopes from —0.07 to —0.04 
dex kpc”*^. The last two works conclude that O and Cl 
gradients are similar, a result which is consistent with a 
lockstep evolution of both elements. However, some recent 
theoretical works indicate that PNe progenitors can mod¬ 
ify their initial O abundance. Some O enrichment can be 
produced in the th ird dredge-up of asymptotic giant branch 
(AGB) stars (e.g. iKarakas fc Lattanzid l2014h but, on the 
other hand, hot bottom burning in t he most massive pro¬ 
genit ors of PNe can destroy it (e.g. IKarakas fc Lattanziol 
l2014h . The efficiency of both processes would depend on 
the mass of the progenitors stars and their metallicity. In 
any case, the observatio nal evidences of such O varia tions 
in PNe are contradictory. Ip^elgado-Inglada et al.l (l2015l i have 
obtained observational evidences of significant O enrichment 
(by ~ 0.3 dex) in PNe with carbon-rich dust, with progeni- 
tor masses in t he ran ge 1.5—3 Mq. But, on the other hand, 
iMagrini et al.l ll2009l i comparing the O and Ne abundances 
of the PNe population of M33, did not find any evidence of 
systematic dredg e-up of O in the meta l licity range of that 
galaxy. Moreover, iMagrini fc GoncalvesI ll2009l 'l reviewed the 
O/H ratios of PNe and Hll regions in dwarf galaxies of the 
Local Group finding possible evidences of O dredge-up only 
at very low metallicities. In contrast, there is not any theo¬ 
retical prediction that the amount of Cl may change along 
the evolutio n of the low and intermediate stars that produce 
PNe. In fact. lDelgado-Inglada et"^ (l2015fl propose Cl/H ra¬ 
tio as a more reliable proxy for metallicity than O/H in PNe 
abundance studies. Therefore, it would be not unexpected 
to obtain somewhat different gradients for O and Cl using 
PNe abundance data. The aim of this paper is to determine 
the Cl and O abundance for a sample of Hll regions with 
the best available spectra and located across the disc of the 
Galaxy to better explore the nucleosynthetic origin of Gl. 

The structure of this paper is as follows. In ^we de¬ 
scribe the set of H ll region spectra used to derive the nebular 
parameters. In lj3]we present the results of our physical con¬ 
ditions and abundances recalculations, as well discuss the 


Table 1. Sample of Galactic Hll regions. 


Hii region 

Rg 

(kpc) 

Data source 

MI 6 

6.34 

Garcfa-R.oias et al. (2006') 

M8 

6.41 

Garci'a-Roias et al. (2007) 

M17 

6.75 

Garci'a-Roias et al. (2007) 

M20 

7.19 

Garcia-Roias et al. (2006) 

NGG 3576 

7.46 

Garcfa-R.oias et al. (2004) 

M42 

8.40 

Esteban et al. (2004) 

NGG 3603 

8.65 

Garcfa-Roias et al. (2006) 

S 311 

10.43 

Garcfa-Roias et al. (2005) 

NGG 2579 

12.40 

Esteban et al. (2013) 


results pertaining to the Galactic Cl radial gradients, the 
local Cl/0 ratio and the empirical ICF{C^'^) we propose. 
Finally, in ^we summarize our main conclusions. 


2 OBSERVATIONAL DATA 


We have used intensity ratios of several selected emis- 
sion lines of s p ectra of Ga lactic Hll regions publis h ed by 
Esteban^_etaL (l2004l . l2013h and iGarcfa-Roias et al.l 


20051 . 20061 . 2 OO 7 II . The spectra were obtaine d with the Ultra¬ 


violet Visual Echelle Spectrograph HJVES. lD’Odorico et al.l 
I 2 OO 0 II at the Very Large Telescope (VLT). Nowadays, this 
dataset is the deepest and most complete collection of spec¬ 
tra of Galactic H ll regions available. The list of objects, 
their Galactocentric distances and the reference of their pub¬ 
lished spectra are presented in T able H] All Galactoce ntric 
distances have been taken from lEsteban et al.l (|2005l l ex¬ 
cept in the case o f NGG 2579. Eor thi s object we use the 
value adopted by Esteban et al] (l2013fl and calculated by 
iGonetti et al.l l|2007l l. The distances are based on kinemati- 
cal or photometrical m ethods and ha ve typical uncertainties 
of less than 1 kpc (e.g. iRussei il l2003h . We have assumed the 
Sun at 8 kpc from the Galactic Centre. 


3 RESULTS AND DISCUSSION 

Although the source papers of the data used in this work 
present calculations of physical conditions - electron tem¬ 
perature, Te, and density, Ue - and ionic and total abun¬ 
dances of some elements, we prefer to re-calculate all the 
relevant quantities in order to have a homogeneous set of 
values calculated with the same methodology and updated 
atomic data. 


3.1 Physical Conditions 

Physical conditions were calculated mak i ng u se of the 
version 1.0.2 of PyNeb (iLuridiana et al.l l2015h . PyNeb 
is an updated and more versatile version of the nebu¬ 
lar package of iraf written in python language. The 
density-sensitive emission line ratios we have used are 
[On] 3726/3729, [Sii] 6717/6731, [Gliii] 5518/5538 and 
[Ariv] 4711/4740, and the temperature-sensitive ones are 
[Nil] 5755/(6548-f6584), [On] (7319-b7330)/(3726-b3729), 
[Sii] (4069-b4076)/(6717-b6731), [Oiii] 4363/(4959+5007), 
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Table 2. Atomic dataset used for collisionally excited lines. 


Transition probabilities 
Ion and energy levels 


Collisional strengths 


N+ 

0 + 

q2+ 

s+ 

S2+ 

Ar2+ 

Ar3+ 

C1+ 

Cl2+ 

Cl^+ 


^oese_^ischer_&_^chiev (200A) 
^oese_^ischer_&_^chiev (200A) 
,Roese_£ischer_&_^chiev (2004 ); ^tore^_&_Zei£gen (2000 
0odobedov^_e0^ (2009) 
0odobedova_e0^ (2009) 

^endo^ (1980 ; JCauftnai^_&_Sugar (1980 
^endoz^_^_Zei££en (19820 
^endoz^_^_Zeig£en (1983) 

^endoza (19^) ; j<aufoiai^_&_Sugar (1980 
^itzsche_et_^^ (1990 

j^aufinai^_^_Sugar (1980 ; 00s_^^artinson (1980 
^endoz^^_Zmg£en (1982b) 


Taval (•2011! 
JGsielius_et_^l^ (2000 
0tore^_et_0j (2014) 
^^^a^_&_^0sminn^ (2010 
(1990 

^0^is_et_0^ (1990 
Ramsbottoii^_&_0eU (1997) 
Taval (2004! 

0utler_^0ei££en (1989) 
^ei0o^ (19^) 
^0a«s_et_ah (1995) 


[Sill] 6312/(9069+9532), and [Arm] 5192/(7136+7751). 
PyNeb permits to easily customize the atomic datasets. In 
our case we have used the atomic data listed in Table [2l 
Temperatures and densities were determined using an itera¬ 
tive procedure. Firstly, we assume an Initial Te of 10,000 
K to determine rie and this last value is used to calcu¬ 
late Te. The uncertainties of each quantity were obtained 
performing Monte Carlo calculations changing the values 
of the emission line ratios within the line intensity errors 
quoted in the source papers of the observational data (see 
Temperatures of low ionization species - re([Oll]), 
re([Nll]) and re([Sll]) - were derived considering the av¬ 
erage of ne([Oll]) and ne([Sll]). For calculating re([Olll]), 
re([Slll]) and re([Arlll]) we considered the average of all 
the density determinations available - ne([Oll]), ne([Sll]), 
ne([Cllll]) and rie ([Ariv]) - because all the density-sensitive 
ratios give similar values and also because ne([Oll]) and 
ne([Sll]) show the lowest uncertainties. Using an average of 
rie ([Cl III]) and ne([Arlv]) would provide almost identical re¬ 
sults but with larger uncertainties. TablesUand [S]show the 
values of electron temperatures and densities, respectively, 
obtained from the different individual line-intensity ratios. 
Considering that Cl^^ is the the dominant ionic species of 
Cl in our Hll region sample and that there are several sets 
of atomic data available, we decided to use combinations 
of those datasets for deriving the rie ([Cl ill]) as well as the 
Cl^’'’/H''’ ratio. In Table [3] we present the different combina¬ 
tions used for this ion. The value of rie ([Cl ill]) included in 
Table [S] is the average of the densities obtained with the six 
combinations. The final adopted electron density for each 
object, <ne>, is the weighted mean of the values obtained 
for the different density-sensitive indicators and it is also 
included in Table [5] 

We assume a two-zone approximation for the nebula 
estimating the representative values of the electron temper¬ 
ature for the zones where low and high-ionization poten¬ 
tial ions are present, re(low) and Te(high). Those values 
have been used to determine ionic abundances, re(low) was 
calculated as the mean of re([Nll]), re([Oll]) and re([Sll]) 
weighted by their uncertainties. On the other hand, re(high) 
was calculated as the weighted mean of re([S m]), Te([Ar iii]) 
and re([Olll]). The corresponding values of re(low) and 
Te(high) used for each object are included in Table [I] 


Table 3. Combinations of atomic datasets used for CU”*". 


Set No. 

Transition probabilities 
and energy levels 

Collisional strengths 

1 

Mendoza ('1988') 

Butler Zeiooen ('1989) 

2 

Mendoza flQSS) 

Mendoza fl983) 

3 

Mendoza flQSS) 
Kaufman & Suerar (1986! 

Butler & ZeioDen ('1989) 

4 

Mendoza fl983) 
Kaufman & Suear ('1986) 

Mendoza fl983) 

5 

Fritzsche et al. fl999') 

Butler & ZeioDen ('1989) 

6 

Pritzsche et al. fl999') 

Mendoza ('1983) 


3.2 Ionic and total abnndances 

Ionic abundances of Cl"^, Cl^"^ and Cl®'*' have 

been derived from CELs under the two-zone scheme using 
the PyNeb package and the atomic dataset indicated in 
Table [2l We have used Te(low) for calculating O’*" and Cl"^ 
and re(high) for O®'*', Cl*^^ and Cl®'*'. Although there are 
several datasets of atomic data for O’*' and O'^'*' available 
in the literature, we have selected those that, according to 
our experience, provide more reliable results. In the case of 
Cl®'*', considering that a) it is the dominant ionic species 
of Cl in our Hll region sample, and b) the availability of 
several references of atomic data, we decided to calculate 
the Cl'^’^/H'^ ratio using different combinations of datasets, 
which are included in Table [3] and indicated with different 
numbers. In contrast, for Cl^ and Cl®'*' we have a single 
option of atomic dataset available and therefore no other 
combinations are possible. 

In Table[6]we present the Cl'^^/H'^ obtained for the dif¬ 
ferent combinations of atomic data used and the weighted 
mean value that we assume as representative for the abun¬ 
dance of this ion for each Hll region. In Table [7] we include 
all the ionic abundances we calculate for the different ionic 
species of Cl and O which emission lines have been detected 
in the nebulae. The uncertainties of ionic abundances were 
obtained performing Monte Carlo calculations. To do that, 
we generate 500 random values for each line intensity using 
a Gaussian distribution centered in the observed line inten¬ 
sity with a sigma equal to the associated uncertainty. For 
higher number of Monte Carlo simulations, the errors in the 
computed quantities remain constant. 
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Table 4. Electron temperatures (K). 


Object 

TedNll]) 

FedOll]) 

Teds II]) 

Te (low) 

Te([OlIl]) 

Te([SlIl]) 

Te([Ariii]) 

Te(high) 

M16 

8360 ± 150 

8330 ± 1210 

9290 ± 730 

8400 ± 190 

7640 ± 140 

8300 ± 270 

— 

7780 ± 170 

M8 

8460 ± 130 

9370 ± 2410 

9600 ± 740 

8520 ± 210 

8060 ± 100 

8460 ±170 

7520 ± 310 

8120 ± 130 

M17 

8890 ± 250 

9000 ± 5130 

8290 ± 650 

8810 ± 310 

8000 ± 120 

8070 ± 230 

8380 ± 410 

8040 ± 160 

M20 

8340 ± 160 

8560 ± 1670 

9050 ± 500 

8410 ± 200 

7870 ± 260 

8080 ± 210 

± 

7990 ± 250 

NGC 3576 

8680 ± 210 

10090 ± 1850 

11280 ± 1620 

8740 ± 250 

8460 ± 130 

8650 ± 400 

8570 ± 360 

8490 ± 180 

M42 

10290 ± 300 

10530 ± 1770 

13610 ± 3550 

10320 ± 360 

8350 ± 90 

9520 ± 580 

8290 ± 260 

8370 ± 120 

NGC 3603 

11360 ± 600 

18980 ± 4370 

18870 ± 3090 

11760 ± 760 

9080 ± 150 

9670 ± 360 

- 

9170 ± 180 

S 311 

9370 ± 210 

10330 ± 1950 

9610 ± 500 

9420 ± 270 

8980 ± 110 

9470 ± 280 

8760 ± 970 

9040 ± 140 

NGC 2579 

- 

12460 ± 1690 

12240 ± 1230 

12340 ± 1280 

9370 ± 60 

10400 ± 290 

8640 ± 290 

9380 ± 80 


Table 5. Electron densities (cm ®). 


Object ne([Sll]) ne([Oll]) ne([Cllll]) rie ([Ariv]) <ne> 


M16 

M8 

M17 

M20 

NGC 3576 
M42 

NGC 3603 
S 311 
NGC 2579 


1100 ± 200 
1280 ± 160 
780 ± 180 
270 ± 80 
1090 ± 320 
4590 ± 3910 
3000 ± 1260 
300 ± 90 
920 ± 170 


1140 ± 230 
1620 ± 630 
500 ± no 
260 ± 50 
1700 ± 380 
6590 ± 1970 
2600 ± 650 
280 ± 80 
1240 ± 280 


1215 ± 560 
1930 ± 320 
250 ± 240 
320 ± 290 
3180 ± 710 
7470 ± 750 
4850 ± 1710 
500 ± 420 
1790 ± 550 


1990: 


2930 ± 1450 
4870 ± 1090 
928: 

1030: 


1120 ± 240 
1420 ± 210 
530 ± 140 
260 ± 60 
1570 ± 410 
6590 ± 1020 
2910 ± 870 
290 ± 90 
1060 ± 220 


The total abundances of Cl and O have been determined 
by simply adding their ionic abnndances and are included in 
the last two columns of Table [T] The four nebulae of lower 
ionization degree of the sample do not show detections of 
[Cl iv] lines and we do not expect any relevant contribution 
of Cl®^ to the total Cl abundance. In Table [8] we compare 
the total Cl abundance we obtain from the sum of Cl^, 
Cl^"*" and Cl^"*" abundances with that we obtain only adding 
Cl^ and Cl^'*' and also the ionization degree of the nebu¬ 
lae (parametrized by the /O^ ratio). As we can see, the 
contribution of the Cl®^/H+ ratio to the total Cl abundance 
is always very small even in the objects of the highest ion¬ 
ization degree and never larger than 0.02 dex (see column 
3 of Table [8]). Therefore, the contribution of Cl®"*" in the 
four Hll regions with the lowest ratios: M20, M16, 

M8 and S 311, should be negligible. For these four objects, 
the assumption of Cl/H = C1'’"/H+-|-CP^/H+ seems to be 
perfectly applicable. 

3.3 The Galactic Cl/H gradient 

We have performed least-squares linear fits to the galac- 
tocentric distance of the Hll regions, Rq, and their corre¬ 
sponding Cl/H, O/H and Cl/O ratios. For the fitting, each 
elemental ratio has been weighted by the inverse of its un¬ 
certainty. These fits give the following radial gradients: 

12 + log(Cl/H) = 5.40(±0.11) - 0.043(±0.012) x Rg\ (1) 
12 -f log(0/H) = 8.83(±0.07) - 0.043(±0.009) x 7 ?g; (2) 

log(Cl/0) = -3.41(±0.12) - 0.001(±0.014) x Rq. (3) 
In Figure [T] we represent the ionized gas phase Cl/H, 


5.4 



4.6 t-1-^-1-1-^^-1-^^-1-1-L 

6 8 10 12 



6 8 10 12 


Galactocentric distance (kpc) 

Figure 1. Ionised gas phase Cl, O and Cl/O abundances of 
the Galactic H II regions studied in this work. The lines indi¬ 
cate the radial gradients obtained from the linear fit of the data. 
The solar symbol ind icates the meteoritic Cl abundance given by 
iLodders et al.l ||2009| '). The Sun is located at 8 kpc. 
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Table 6. Cl^+/ H+ ratios® determined from combinations of atomic data sets and adopted mean value 


Object 


1 


2 


3 


4 


5 


6 

Adopted 

MI6 

5.06 

± 

0.04 

5.08 ± 0.04 

5.05 

± 

0.04 

5.08 ± 0.04 

5.06 

± 

0.04 

5.08 ± 0.04 

5.07 

± 0.04 

M8 

5.03 

± 

0.02 

5.05 ± 0.02 

5.02 

± 

0.02 

5.04 ± 0.04 

5.02 

± 

0.02 

5.05 ± 0.03 

5.04 

± 0.03 

M17 

5.05 

± 

0.03 

5.07 ± 0.03 

5.05 

± 

0.03 

5.07 ± 0.03 

5.05 

± 

0.03 

5.07 ± 0.03 

5.06 

± 0.03 

M20 

4.99 

± 

0.04 

5.01 ± 0.04 

4.99 

± 

0.04 

5.01 ± 0.04 

4.99 

± 

0.04 

5.02 ± 0.04 

5.00 

± 0.04 

NGC 3576 

4.93 

± 

0.03 

4.96 ± 0.04 

4.93 

± 

0.04 

4.95 ± 0.05 

4.93 

± 

0.04 

4.96 ± 0.05 

4.94 

± 0.04 

M42 

5.12 

± 

0.03 

5.14 ± 0.02 

5.11 

± 

0.03 

5.13 ± 0.05 

5.11 

± 

0.02 

5.13 ± 0.03 

5.12 

± 0.03 

NGG 3603 

5.02 

± 

0.05 

5.05 ± 0.05 

5.02 

± 

0.06 

5.04 ± 0.06 

5.02 

± 

0.05 

5.04 ± 0.06 

5.03 

± 0.05 

S 311 

4.85 

± 

0.03 

4.88 ± 0.03 

4.85 

± 

0.03 

4.87 ± 0.03 

4.85 

± 

0.04 

4.88 ± 0.03 

4.86 

± 0.03 

NGG 2579 

4.80 

± 

0.03 

4.83 ± 0.03 

4.80 

± 

0.03 

4.82 ± 0.04 

4.80 

± 

0.03 

4.83 ± 0.03 

4.81 

± 0.03 


® In units of 12+log(Cl2+/H+). 


Table 7. Ionic and total abundances for O and Cl® 


Object 

0 + 

02+ 

G1+ 

CI 2 + 

gp+ 


0 



G1 


MI 6 

8.44 

± 

0.04 

7.88 

± 

0.03 

4.42 

± 

0.04 

5.07 

± 

0.04 


— 


8.54 

± 

0.04 

5.15 

± 

0.04 

M8 

8.32 

± 

0.06 

7.87 

± 

0.02 

3.99 

± 

0.23 

5.04 

± 

0.03 


- 


8.45 

± 

0.04 

5.09 

± 

0.03 

M17 

7.81 

± 

0.05 

8.43 

± 

0.03 

3.59 

± 

0.10 

5.06 

± 

0.03 

3.10 

± 

0.22 

8.53 

± 

0.02 

5.08 

± 

0.03 

M20 

8.44 

± 

0.04 

7.68 

± 

0.04 

4.44 

± 

0.04 

5.00 

± 

0.04 


- 


8.51 

± 

0.04 

5.11 

± 

0.04 

NGG 3576 

8.09 

± 

0.06 

8.36 

± 

0.03 

3.80 

± 

0.06 

4.94 

± 

0.04 

3.21 

± 

0.05 

8.55 

± 

0.03 

4.98 

± 

0.04 

M42 

7.75 

± 

0.06 

8.42 

± 

0.02 

3.49 

± 

0.07 

5.12 

± 

0.03 

3.65 

± 

0.05 

8.50 

± 

0.02 

5.15 

± 

0.03 

NGG 3603 

7.29 

± 

0.09 

8.41 

± 

0.03 

3.13 

± 

0.06 

5.03 

± 

0.05 

3.85 

± 

0.03 

8.44 

± 

0.03 

5.06 

± 

0.05 

S 311 

8.25 

± 

0.05 

7.81 

± 

0.02 

4.29 

± 

0.04 

4.86 

± 

0.03 


- 


8.39 

± 

0.03 

4.96 

± 

0.03 

NGG 2579 

7.35 

± 

0.12 

8.20 

± 

0.02 

3.25 

± 

0.07 

4.81 

± 

0.03 

2.78 

± 

0.05 

8.26 

± 

0.03 

4.83 

± 

0.03 


® In units of 12+log(X'+/H+) or 12+log(X/H). 


Table 8. Contribution of Cl®"*'/H+ to the Cl/H ratio. 


Object 

Cl® 

Cl+l-CP+b 

dif.‘= 

log{0^+/0+) 

M20 

5.11 

5.11 

— 

-0.76 

M16 

5.15 

5.15 

- 

-0.56 

M8 

5.09 

5.09 

- 

-0.45 

S 311 

4.96 

4.96 

- 

-0.44 

NGG 3576 

4.98 

4.97 

0.01 

0.27 

M17 

5.08 

5.07 

0.01 

0.62 

M42 

5.15 

5.13 

0.02 

0.67 

NGG 2579 

4.83 

4.82 

0.01 

0.85 

NGG 3603 

5.06 

5.04 

0.02 

1.12 


® In units of 12+log(Cl/H). 
b In units of 12+log(Cl+/H++Cl2+/H+). 
Difference between columns 2 and 3. 


Table 9. Determinations of the Galactic Cl/H abundance gradient. 


0/H and Cl/0 ratios of the Hll regions and the corre¬ 
sponding radial abundance gradients given above. In Table[9] 
we compare our Cl gradient with previous ones determined 
from PNe data. The slope of our Cl/H gradient is somewhat 
flatter tha n those determined by Faundez-Abans fc Maciell 
lIlQBitl and iMaciel fc Chiapninil (|l994h . This may indicate an 
evolutionary effect, because a temporal flatting of the radial 
abundance gradients h as been predicted in different chemica l 
evolution models (e.g. IHou et aDl200fll : iMolla fc DfadfioObl l . 
H II regions abundances represent the present-day composi¬ 
tion of the ISM but, in PNe, abundances of most elements 
should reflect the ISM composition in the past, at the mo¬ 
ment of the formation of their progenitor stars, which have 
a wide age span, from about 1 Gyr to about 8 Gyr (e.g. 
iMaciel et al.ll2009l 'l. However, this possible temporal flatten¬ 
ing of the Cl gradient is doubtful because the parameters of 
our fit shown in Figure [T]- slope and in tercept - are in ver y 
good agreement with those obtained bv iHenrv et al.l ll2004h . 


Reference 

Objects 

Slope 

(dex kpc”^) 

Intercept 

F-AM87® 

PNe 

-0.057 

5.67 


PNe 

-0.07 ± 0.01 

— 

HKB04‘= 

PNe 

-0.045 ± 0.013 

5.45 ± 0.11 

This work 

Hii regions 

-0.043 ± 0.012 

5.40 ± 0.11 


® iFaundez-Abans &: MacidI lll987f) . 
^ iMaciel fc Chiappinil |l99j. 
iHenrv et al.l ll2004h . 


As we can see comparing the slopes of the gradients 
we obtain, the Cl/H and O/H gradients are identical within 
the uncertainties and this is reflected in a constant Cl/0 
ratio across the galactic disc. This is consistent with a lock- 
step evolution of both elements. As it w as commented in m 
a sim ilar c onclusion wa s obta ined by IMaciel fc Chiappinil 
lll994h and Henry_et_alJ (12004) from large datasets of PNe 
spectra. IClavtonI (I 2 OO 3 I) and lidenry et all (l2004l ) comment 
that theoretically a certain amount of Cl may be produced 
by Type I supernova, but our results indicate that this 
amount should be negligible compared to the contribution 
due to massive stars. Theoretical predictions also indicate 
that O should also be in part produced by Type I super- 
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nova, bu t their relative cont ribution is only of the order of 
8% fe.g. IWooslev et al]|2002l '). 

In Figure [T] we can see that two objects show somewhat 
discrepant Cl/H and Cl/O ratios with respect to the others. 
One is NGC 3576, at Rg = 7.46 kpc, and the other M42, the 
Orion nebula, at Rg = 8.40 kpc. In the case of NGC 3576, 
the rather low Cl/O ratio is produced by the combination of 
a low Cl/H and a somewhat high value of O/H. In the case 
of M42 the reason seems to be basically a higher value of 
Cl/H. It is difficult to ascertain the reason of such differences 
but they might be due to relative different contribution of 
local sources of O and Cl in the chemical composition of the 
nebulae or a different relative depletion fraction onto dust 
particles in those two objects. 


3.4 Comparison with the solar Cl abundance and 
the local Cl/O ratio 


lAsolund et al.l l|2009l l indicate that the solar photospheric Cl 
abundance is very uncertain, 12+log(Cl/H) = 5.50 ± 0.30. 
This value cannot be derived from spectroscopy of the quiet 
Sun but, instead it was de termined from IR spectrum of 
sunspots using lines of HCl ( Hall_&No^ Il972ll . Meteoritic 
values presented by iLodders et al.l ( 2009 ) are more confi¬ 
dent, they give a value of 5.25 ± 0.06, which is included in 
Figure[T]for comparison. From our radial Cl/H gradient, we 
estimate a ionized gas phase 12 -|- log(Cl/H) = 5.06 at the 
solar Galactocentric distance of 8 kpc. Therefore, the mete¬ 
oritic Cl abundance is of the order of 0.2 dex higher than 
the nebular one. This is an apparently puzzling situation 
because meteorites should represent the composition of the 
solar nebula ~4.5 Gyr ago and one would expect that the 
Cl abundance at the present-day ISM to be higher than the 
meteoritic value, just the opposite of what we find. Since 
- as far as we know - there are not numerical determina¬ 
tions about the evolution of Cl in the disc of the Milky Way, 
one would expect some increase of the Cl/H ratio in the 
last ~ 5 Gy r. In fact, chem i cal ev olution models by ICarigil 
l|2003h and lAkerman et al.l l|2004h estimate that the O/H 
ratio at the solar Galactocentric distance has increased by 
0.12 dex since the Sun was formed. It is important to remark 
that the same problem exists when one compares the solar 
O/H ratio with that of ionized nebulae in the solar neigh¬ 
borhood. The solar photospheric O/H ratio is about 0.2 dex 
higher than the nebular values when determined from the in¬ 
tensity of collisionally excited lines (CELs) of ionic species 
of O. The reason of such difference can be due to deple¬ 
tion onto dust grains in the Hll regions and/or physical 
mechanisms associ ated to the so-called abundanc e discrep¬ 
ancy problem fe.g. Gin:cfaj:Rojafi__&J] steb.Mj| 20 ^ Accord- 


ancy problem le.g. Hja rcia-Ko ias Jj^s tebanllZ UUill. Accorg- 

ing tolMesa-Delgado et al.l(2009lj andiPeimbert fc Peimberd 


ll2010ll . the fraction of O embedded in dust grains in the 
Orion nebula is about 0.12 dex and, therefore, can only ac¬ 
count for part of the problem. On the other hand, the abun¬ 
dance trend of Cl in diffuse and dense interstellar clouds 
is rather irregular and amount of dust deple tion can not 
be estimated for this element lljenkind l2009l l. The abun¬ 
dance discrepancy problem arises when ones compares abun¬ 
dances of a given ion obtained from CELs and recombina¬ 
tion lines (RLs). In all the Hll regions included in this pa¬ 
per, O^^ abundances determined from CELs are between 
0.2 and 0.3 dex higher than those determined from RLs 



0 / 0 ^ + 


Figure 2. Values of as a function of the 0/0^"^ ratio 

for our H II regions sample. The continuous line represents a linear 
fit to the data. 


(see ICarcfa-Roias fc EstebanI l2007l l . Among different pos¬ 
sibilities, one of the mechanisms that may produce such dis¬ 
crepancy is the presence of fluctuations in the spatial dis¬ 
tribution of electron temperat ure in the nebu lae (the so- 
called temperature fluctuations. |Peimbertl[T967^ . Under the 
presence of temperature fluctuations, ionic abundances de¬ 
termined from CELs would be lower than the real ones. A 
qualitatively similar effect would produce the presence of 
a “kappa-distribution” in the energy spectrum of free elec¬ 
trons in Hll regions, a nother explanation f or the abundance 
discrepancy problem llNicholls et al.l[2012ll . In the presence 
of any of both mechanisms, temperature fluctuations or a 
“kappa-distribution”, the effect on the derived Cl/H would 
be similar, the Cl/H ratio we have determined in this pa¬ 
per - from the intensity of CELs - would be lower than the 
real ones, and therefore it can explain the difference between 
solar and nebular Cl abundances. 

From our data, we have determined the weighted mean 
of the Cl/O ratio across t he Galactic disc a nd this is 
log(Cl/0) = —3.42 ± 0.06. iHenrv et al.l (l2004l ~) compile a 
number of determinations of the Cl/O ratio based on PNe 
and Hll region observations, and the logarithmic values 
range from —3.67 to —3.30. Therefore, our mean ratio is 
consistent with those previous calculations. Moreover, from 
the list of recommended present-day solar ab undances from 

S spheric and meteoritic data compiled bv iLodders et al] 
, we obtain a log(Cl/O )0 = —3.48 ± 0.09. In rather 
good agreement with our nebular determination. This simi¬ 
larity and the rather small dispersion of the individual val¬ 
ues of the Cl/O ratio we obtain for the different Hll regions 
suggest that the individual or combined action of the two 
aforementioned effects - dust depletion and abundance dis¬ 
crepancy - should affect similarly to both elements - Cl and 
O - in all objects. 


3.5 An empirical ionization factor for Cl 

Usually, [Cl III] lines are the only emission lines of Cl that 
can be measured in the spectra of Hll regions. Therefore, 
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Table 10. Comparion of Cl abundances determined with different ICF schemes. 


Object 

0/02+ 

This work 

12 -t log(Cl/H) 

P&T-P77=>- M&R91*" ISMGT06“ 

D-IMS14'l 

M20 

6.76 

5.11 

5.18 

5.12 

5.18 

5.17 

M16 

4.57 

5.15 

5.33 

5.24 

5.23 

5.22 

M8 

3.80 

5.09 

5.23 

5.13 

5.18 

5.18 

S 311 

3.80 

4.96 

5.12 

5.01 

5.00 

5.00 

NGC 3576 

1.55 

4.98 

- 

5.18 

4.99 

5.06 

M17 

1.26 

5.08 

- 

5.42 

5.14 

5.22 

M42 

1.20 

5.15 

- 

5.54 

5.22 

5.29 

NGC 2579 

1.15 

4.83 

- 

5.52 

4.98 

5.01 

NGC 3603 

1.07 

5.06 

- 

5.52 

5.36 

5.28 


^ IPeimbert k. Torres-Peimberd |l973). 
^ iMathis &: Rosal lll99ll V 
^ I Izotov et al.l ll2006h . 
iDelgado-Inglada et al.l ll20l4l 


an ionization correction factor (ICF) has to be adopted to 
estimate the total Cl abundance. The ICF{C\^'^) is a multi¬ 
plicative factor to transform Cl^^/H^ ratios into Cl/H ones 
following the simple relation: 

Cl/H = ICF{Cf+) X Cl^+/H+. (4) 

jPeimbert fc Torres-Peimberd (Il977h proposed a relation 
based on the similarity of the ionization potential of the 
ionic species of Cl with those of O and S. However, they 
only re commend its u s e for values of ICF{C\-^^) close 
to one. iMathis fc Rosal |l99j) presented a quite complete 
ICF scheme for several elements based on nebular pho- 
toionlzation models using mainly plane-parallel stellar at¬ 
mospheres as ionization sources. These authors admit that 
their ICF{C\-^^) were very uncertain, with abundance dif¬ 
ferences of factor of two f or a given obj e ct dep ending on the 
model atmosphere used, lizotov et al.l (l200m obtained fit¬ 
ting expressions based on photoionization models for extra- 
galactic H II regions and defined diffe rent fittings depending 
on th ree metallicity ranges. Finally, iDelgado-Inglada et al. 
ll2014li compute a large grid of photoionization models cov¬ 
ering a wide range of physical parameters representative of 
PNe. This work includes an ICF( Cl^'^/O^) that can be 
appli cable for H ll region conditions (IDelgado-Inglada et al.l 
l2015l i. In Table fTOl we compare our direct determinations of 
the Cl/H ratios for the Hll regions sample with the val¬ 
ues we obtain using our quoted Cl^'''/H''' ratios and the 
ICF schemes mentioned above. In Table [IQ] we also in¬ 
clude the O/O^”*" ratio as indicator of the ionization de¬ 
gree of the objects. In the case of the determination s base d 
on the ICF scheme by IPeimbert fc Torres-PeimbertI lll977ll . 
we have only included objects with ICF< 2.0. In 
the case of th e abu ndances obtained using the ICF of 
IMathis fc Rosal (1 199 ill . we use the coefficients for “cool” Ku- 
rucz atmospheres for all the objects except for NGC 2579 
and NGC 3603, for which “hot” Kurucz atmospheres are 
more appropriate. The inspection of the values included in 
Table fTUl indicates that all the ICF schemes provide Cl/H 
ratios systematically larger than our total Cl abundance 
determinations. All the schemes overestimate the contribu¬ 
tion of Cl*^ at lower ionizatio n degrees but also of Cl^"*" at 
higher ones. The correction by IPeimbert fc Torres-PeimbertI 
lll977h . even of a very restricted applicability, gives values 


between 0.1 and 0.2 dex higher than our direct determi¬ 
nations. This is the /C'F’(CP^) that gives the larger dis¬ 
crepancies for the low ionization objects. The rest of the 
ICF schemes give differences smaller than 0.1 dex for the 
low ionization nebulae, those having 0/0^"^ > 2. The situ¬ 
ation becomes worst in the case of the high ionization Hll 
regions, the objects with 0/0 ^"*" < 2. In this regime, the 
ICF of IMathis fc Ros^ lll991^ give differen ces larger than 
0.2 dex. The scheme of lizotov et al.l (l2006fl gives the best 
results, with differences smaller that 0.2 dex - even smaller 
than 0.1 dex for most of the objects - except in the case of 
N GC 3603, the object wit h the highest ionization. The ICF 
of IDelgado-Inglada et al] ll2014ll also gives reasonably good 
results with differences of the order or lower than 0.2 dex 
in all the cases. It is clear that all the ICF schemes over¬ 
estimate the contribution of Cl®'*' and that this difference 
increases with the ionization degree of the nebulae. 

Our results lead us to propose a new empirical 
ICF{Cl^'^) based on observations and valid for the ioniza¬ 
tion range covered by our objects. In Figure[2]we present the 
value of the ICF{Cl^'^) we obtain for each object - which 
is simply the ratio between the total Cl abundance with re¬ 
spect to the Cl®'*' abundance - versus the ionization degree 
of the nebulae, parametrized by the 0/0®’*' ratio. The hg- 
ure indicates a rather clear correlation that suggests nothing 
more complicated than a simple linear relation. The least- 
squares linear fitting is: 

ICF{Cf+) = 1.02 -f 0.04 X 0/0®-^. (5) 

The uncertainty associated with our ICF is about 
0.03 dex. The range of validity of this ICF{Cl^'^) correspond 
to the interval of O/O®"*" ratio covered by our objects: 1 < 
0/0®'*' < 7. 


4 CONCLUSIONS 

We determine the radial abundance gradient of Cl of the 
Milky Way from published spectra of H ll regions. The 
dataset is the deepest and most complete collection of spec¬ 
tra of Galactic H ll regions available. We have re-calculated 
the physical conditions - electron temperature and density - 
and ionic and total abundances of G1 and O using the same 
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methodology and updated atomic data for all the objects. 
The total abundances of Cl and O are determined by simply 
adding their ionic abundances, without assuming ionization 
correction factors. 

We find that the slopes of the radial gradients of Cl 
and O are identical within the uncertainties: —0.043 dex 
kpc“^. This is consistent with a lockstep evolution of both 
elements. There are previous determinations of the Galac¬ 
tic gradient of Cl but from observations of PNe. However, 
this is the first time that the gradient is computed from 
direct - not ICF based - abundance determinations of 
Cl in Hll regions. We obtain that the mean value of the 
Cl/O ratio across the Galactic disc is log(Gl/0) = —3.42 ± 
0.06, consi stent with previous determinations based on PNe 
data (e.g. iHenrv et al.l l2004ll and from the recommended 
solar photospheric a nd meteoritic abundances compiled by 
iLodders et al.l ll2009h . 

We compare our direct determinations of the Cl/H ra¬ 
tios with those determined from our Gl^^ abundances and 
using some available ICF schemes of the literature. We find 
that our total G1 abundances are always lower than the val¬ 
ues determined using ICFs, indicating that those correc¬ 
tions schemes systematically overestimate the contribution 
of G1+ and Gl®+ species to the total G1 abundance. Finally, 
we propose an empirical ICF{Cl^'^) to estimate the total 
Cl abundance in Hll regions when only the Cl^"*"/ ratio 
can be determined from their spectra. 
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